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We show that the CJT thermodynamic potential of dynamical quark models with a quark propa- 
gator represented by complex conjugate mass poles inevitably exhibits thermodynamic instabilities. 
We find that the minimal coupling of the quark sector to a Polyakov loop potential can strongly 
suppress but not completely remove such instabilities. This general effect is explicitly demonstrated 
in the framework of a covariant, chirally symmetric, effective quark model. 
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Q-i I. INTRODUCTION 

The fate of hadronic matter in extreme environments, e.g., in the interior of compact stars or 
in the early universe, remains one of the most interesting unanswered questions today. Heavy-ion 
experiments performed at the Relativistic Heavy-Ion Collider (RHIC) at the Brookhaven National 
Laboratory or at the Large Hadron Collider (LHC) at CERN Geneva show that at sufficiently 
high temperature or density, hadronic matter undergoes a phase transition by dissolving into its 
constituents: quarks and gluons. Details of this transition are encoded in the phase diagram of 
quantum chromodynamics (QCD). 

Due to its non-perturbative nature at low temperatures and densities, QCD is best studied on 
the lattice in this regime. Results for the QCD equation of state (EoS) have recently become 
available also at physical quark masses and were extrapolated to the continuum [IHH]. However, 
at finite densities lattice simulations are still limited due to the known sign problem. 

In this situation, effective models of QCD serve to interpret and also extrapolate lattice re- 
sults. To be realistic, such a model must dynamically break chiral symmetry and confine the 
colored degrees of freedom. In particular, quark and gluon propagators should strongly differ from 
their high-energy counterparts. A suitable continuum, non-perturbative and covariant approach is 
provided by the Dyson- Schwinger equations (DSE) (for reviews see e. g. [1-Qj]]) and its descen- 
dants, the non-local chiral quark models fl2l - [T5j (see also (l6l - [l8| b most recently augmented by 



the Polyakov loop (PL) [19|-|2e 

A non-perturbative quark propagator is a solution of its DSE, within the appropriate symmetry- 
preserving truncation scheme. The most crucial behavior of these continuum studies is a strong 
infrared running of the quark dynamical mass, which is to be interpreted as dynamical chiral 
symmetry breaking (DChSB), and of the wave function renormalization. Confinement can be 
realized through the absence of poles of the quark propagator at real time-like four-momenta by 
the criteria of positivity violation 0, Indeed, it is a common feature that, due to strong 

dynamics, the poles are pushed deep into the complex four-momentum region [27H32j |. The most 
simple resulting structure is a series of quartets of complex conjugate mass poles (CCMPs). 
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For a successful phenomenological study it is sufficient to model the strong interactions via an 
effective gluon propagator to be fixed by the infrared observables in the vacuum. Enhancement of 
the interaction in the infrared pushes the quark poles away from the real axis. Another possible 
realization is the absence of propagator poles in the entire complex plane which can be realized 
either by the presence of cuts instead of poles 1331 or by the absence of both when the quark 
propagator is obtained as an entire function (3~3. |35|. A non-constant entire function must have 
a singularity at infinity. In order for such a quark state to become deconfined and to restore its 
approximate chiral symmetry, e.g., at high temperature, the gluon sector of the theory must be 
restructured in the transition region in order to allow for the appropriate changes in the analytic 
properties of the propagators such as the appearance of quasiparticle poles. But in that case, 
the usual strategy is no longer applicable: to predict the behavior at finite temperature and 
chemical potential from a straightforward generalization using the Matsubara formalism without 
changing the analytic properties which were adjusted by constraints from vacuum observables. A 
modification of this strategy lies beyond the scope of the present study. 

A CCMP structure in the quark propagator is sufficient to ensure violation of reflection positivity 
[3(1 |H, HI] , and as such provides a useful form to fit the lattice quark propagator [H, [36| . This 
led to applications at finite quark chemical potential [37j , or for parton distribution functions [38j . 
On the other hand, QCD bound states are affected by this structure [Hj], e.tr a sufficiently heavy 
meson state (typically of the order of I GeV) has unphysical dq thresholds [4(| if a naive analytic 
continuation to the mass pole of the bound state is used (see [12j for a suggestion how this problem 
could be circumvented by a more elaborate analytic continuation). 

In this work we want to further investigate the properties of models with CCMPs by concentrat- 
ing on the finite temperature and finite quark chemical potential EoS. We find that the pressure 
of quark matter leads to an unsatisfactory scenario; the EoS exhibits oscillations in temperature 
which are in turn a consequence of the imaginary part of the CCMPs. We provide analytic insight 
into the nature of these oscillations, and suggest a partial solution to this problem by coupling the 
system to the Polyakov loop (PL). 

This paper is organized as follows. In Section II, we present our arguments in a general form by 
postulating a CCMP parametrization of the quark propagator for which we then obtain the kinetic 
part of the QCD partition function in the quark sector. This result elucidates that the presence 
of complex conjugated mass poles in the quark propagator entails thermodynamic instabilities. 
We introduce the PL variable in the partition function and show that this step very effectively 
suppresses the instability. A separate analysis is performed for the EoS at zero temperature where 
for the CCMP parametrization the quark number density and the pressure can be obtained in closed 
form. In Section III we present the example of dynamical quark models with chirally invariant 
nonlocal interaction, including explicit numerical results and their discussion. In Section IV we 
give the summary and conclusions of our study. 

II. THERMODYNAMICS IN THE CCMP REPRESENTATION 

The (unrenormalized) thermodynamic potential for the quark sector of QCD can be given in the 
form of the Cornwall- Jackiw-Tomboulis (CJT) effective action (4fJ as 

n(T )fi ) = T[S] = -TrLog(5- 1 ) + Tr[ES] + *[S] , (I) 

where S" -1 represents the inverse of the full quark propagator in Euclidean space, 

S-'iPn) = i(7 ■ P)A(Pl) + i-H&nCffi) + B(pl) . (2) 

The quark dressing functions A{p^), B(pf t ) and C(p^) encode effects of the quark selfenergy 
£ = S^ 1 — Sq 1 , expressing all deviations from the free propagator Sq 1 due to nonperturbative 
interaction effects. At finite temperature and chemical potential the "shifted" fermionic Matsub- 
ara frequencies w„ = uj n — ifi = (2n + l)nT — i/j, with temperature T and chemical potential /it, 
are introduced, so that p n = (p, (!>„)• For the free quark propagator Sq we have A = C = 1 and 
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B = m, the current quark mass. The Tr operation implies summation over internal degrees of 
freedom and p n . The leading term in the 2PI loop expansion is the usual one-loop contribution 
which we denote as the kinetic contribution to the thermodynamic potential of the system 

O kin = TrLog(S) , (3) 

and the functional ^[S] contains all higher loop diagrams. In the widely used rainbow-ladder 
approximation (see e.g. Q) it is given by 

= ~Tr\ES\ . (4) 

A quite general and mathematically simple realization of the analytic structure of the quark 
propagator exhibits a series of CCMPs in the p 2 plane for p 2 = —m\ and, p 2 = —ml 2 . The 
complex numbers are ordered such that \mk+i/mk\ > 1. Let us for definiteness assume that 
all the poles are simple. The following arguments can easily be generalized if the quark propagator 
has higher order poles, or branch cuts. At zero three-momentum p = the pole structure in the 
p° = ip4 plane is a series of quartets of poles located at im^ = ztm^ztiml and ±?n£ = ±m^=F«m^., 
and we define rrif. such that m^, mj. > for all k. With p ^ the poles are given by [TJI 

E 2 = p 2 +m 2 . (5) 

where for each k their locations form a quartet in the complex energy plane at i£k = i e fc ± ijk 
and ±££ = ±e fc T *7fc, with 

(ml) 2 + P 2 + VW) 2 - K) 2 + P 2 ] 2 + 4«-) 2 K)4 V2 

J (6) 



The analytic structure of the quark propagator governs the thermodynamical properties of the 
system. Here we perform a simple calculation of the kinetic contribution to the thermodynamic 
potential in the quark sector, with the proposed form of the quark propagator. 



«) 2 - 
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A. Consequences for the quark sector at finite temperature 

Performing the trace in Dirac, color, flavor and momentum space, the kinetic term can be written 



as 

+00 



SlKn(T,») = -2N c N f T [ ^ log [p 2 A 2 (p 2 n )+Q 2 n C 2 (p 2 n ) + B 2 (p 2 n )] . 

For simplicity we work with Nf equal flavors. N c = 3 is the number of colors. 

In order to perform the Matsubara sum we introduce generalized occupation numbers 



(7) 



n± 



(z) = (l + e^ 2 ^)- 1 (8) 



having simple poles at z = iu> n - With the help of the residue theorem, the Matsubara sum is 
converted to an integral along straight lines Re(z) = fi — S and Re(z) = (i + 5, where S > is 
infinitesimal 



dzn + {z)\ogV{z) + / dzn + {z)\ogV{z) 

-ioo~\-fi-\-5 J +'ioo+/j, — 5 

+ oo ^ (9) 

2ni (-T)log[D(iQ n )} 
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where we defined 



V(z) = p 2 A 2 {p 2 ,-z 2 ) - z 2 C 2 {p 2 -z 2 ) + B 2 (p 2 ,-z 2 ) . 



(10) 



Here and in the following we suppress the p 2 dependence in T> for brevity. 

Due to the known analytic structure of the quark propagator the integrals in ([9|) can be evaluated. 
We now close the contour running from — ioo + fj, + 8 to +ioo + fi + S by a large semicircle on the 
positive real axis, and denote this as C\. Then we can rewrite 



V(z) 

T{-2m) {log l + e-^ £k -^ 



I J dzlog U+e-^-^ 



k,ek>fJ. 



+ log 



(ii) 



where the first equality follows from partial integration. The last line is the result of the residue 
theorem, and £k are the previously defined poles. 

For the second term we first make use of the clockwise oriented contour Ci defined as a rectangle 
having vertices in (+ioo + jj, — 5, —ioo + /i — 5, —ioo, +ioo) to obtain 



h=T <b dz log 
'c 2 



1 + e-K'-A 



V{z) 



+ ZOO 



dzn+(z) \ogT>{z) . 



(12) 



If the second term in (fT2|) is rewritten using n + (z) — I — n_(— z), the first of the two resulting 
terms can be Wick rotated to the real axis, providing the vacuum contribution, while the second 
term can be evaluated by yet another contour, defined as C3, where we close the line running from 
—ioo to +ioo by a large semicircle on the negative real axis. This gives 



1 



h = T f dz log 

dp4logT>(ip 4 ) 



Viz) 
V(z) 



T f dz log 
ic 3 



1 + e 



V{z) 



T(-2wi) J2 { lo S [l + e^* - " 3 ] + log [l + e^ 5 *"^] J 



(13) 



k,0<ek<H 



r r 1 r "n f 

+ T(-2Tri) Y { lo S l + e- fj{£k+ ^ +log l + e - /3( ^+^ \-i I dp 4 \ogV(ip A ) 

Collecting the obtained formulas, we can state the kinetic contribution to the pressure 



fi kin (r, ju) = fi zpt - 2TN C N S J2 



k=l 



d 3 p 
(2^)3 



{log 



loe 



I _|_ e -/3(£fc-M) 
1 + e -{S(£ k +v) 



+ log 
+ log 



1 + e-m+f)'] } , 



(14) 



where il zp t is the (divergent) zero-point energy contribution 

d 4 p 



n zpt = -2N c N f 



(2tt)< 



\og[p 2 A 2 (p 2 )+B 2 (p 2 )] 



(15) 



It is plain to see that in the special case of just one pair of real poles at ±m, the original dispersion 
([6|) is reduced to the one of a free relativistic particle, and, accordingly, the second term in Eq. 
(fT4"|) is reduced to the free Fermi gas expression, a situation also encountered, e.g., in the Nambu 
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Jona-Lasinio (NJL) model [42J,|43(, see, e.g., [4J-|47(. In the first term, the integral over p4 can be 
evaluated as well, leading to 



fizpt = -*N c N f J2 



d 3 P (£ £* 
2+T 



iJ (27T) 

which, again in the case of a pair of real poles, is just the usual zero-point energy. 
By combining the logarithms, (I14[) can be cast in a more transparent form 



Q^ in {T,/ J L) = Q zpt -2TN c N f 



log [l + 2e-^ k -^ cos(/3 7fc ) + e ~ mck -^ 



(16) 



+ log [l + 2e-^x + ri cos(/37 fc ) + e -V{^+n)\ } . 



Note that the oscillating cosine functions in the thermodynamic potential could render the quark 
matter unstable. Their origin is traced back to the appearance of imaginary parts "fk of the quark 
mass poles. 



B. Introducing the Polyakov loop 



The traced PL j>(x, T) and its conjugate 3>(x, T) are order parameters for confinement in 
quenched QCD [H, l49j . and as such represent important configurations of the gluon field that 
should be accounted for in the effective thermodynamic description of QCD. They are given as 
thermal expectation values 



(17) 



where V is the untraced PL. In the Polyakov gauge [48j the latter takes a simple form V = 
e*( A 303+A 8 08) ) w here A3 } 8 are color Gell-Mann matrices, with 03 t s being the background gluon field. 
Quark (antiquark) confinement is then signalled by $ = = 0). 

Coupling of the PL to the quarks amounts to a modification of the quark occupation number 
function 



( z ) jl + e /^T(M-i(A 3 <fe+A 8 08))] I 1 



(18) 



Following the same steps as in the previous subsection, the kinetic contribution to the thermody- 
namic potential can be written as 



n kin (T,fi) = -2N f Tj2j ^tr c {log 



log 



loe 



1 + 7?t e -0(**+iO + io g 1 + pt e -«£!+/*) J 



(19) 



Working out the color trace gives 

d 3 p 



n kin (T,/i)=-2iV / T 



log 



(20) 



log [l + We-Wi-ri + 3^e~ 2 ^-^ + e- 3 ^-^] + (jt -> -/x)} . 



Again, in the special case of just one pair of real poles at ±m, this expression is reduced to the 
corresponding term of the Polyakov-Nambu- Jona-Lasinio (PNJL) model (see, e.g., (50l - [53j ). 



G 



The logarithms can be combined to obtain 

d 3 p 



loe 



+ 6$ L-ft**-ti cos(/3 7fc ) + e-^*-^ cos(2/3 7fe )) 



2NfT ^ J (2,)3 

k=l J v ; x 

+6<l L-W(z»-») cos(2/3 7fe ) + e- 8 ^ 6 *-") cos(/3 7fe )) + 9$ 2 e" 2/3(£fc -^ + c^V 4 ^"^ 



+ 18$l'e- 2/3(efc -' i) cos(/3 7fc ) + 2e- 3/3 ( £fc -^) cos(3^ 7fc ) + e -o/9(«s*-/*)l + ( M _> _ M ) 



(21) 



Comparing this with Eq. (|16p. we see that now the dominant cosine terms are weighted by the 
PL. As a consequence, the pressure instabilities are highly suppressed in the confined phase: As 
long as $ and $ are zero, there remains only one cosine term, which is, however, suppressed by 
the third power of the Boltzmann factor. In fact, the mechanism is basically the same as in the 
PNJL model, where the coupling to the PL suppresses the quark degrees of freedom at low T, but 
does not eliminate them entirely [5iT[53l | . 



C. Zero temperature, finite chemical potential 

In this part, special attention is devoted to the effects of the CCMPs along the T = 0, fi > 
axis. The Matsubara sum in ([7]) gets converted to an integral over p 4 . 



fikin(0,/x) = -2N c N f 



d 4 P 



I log[p 2 ^ 2 (p 2 )+p 2 C 2 (p 2 ) + S 2 (p 2 )] 



(22) 



where 



f = v 2 



Pi, P4=P4-1^ 

We start by considering the quark number density 



n{fi) 



2NfN c 



d 4 p ,dV(ip 4 ) 1 
? (-2«p 4 ) 



(27T) 



dpi V{ip4 



(23) 



where T>{ip^) is given by (1101) . At zero chemical potential, the quark density is zero, as it is 
obvious from the integrand being an odd function of p 4 . This allows for the evaluation of the 
integral by a clockwise oriented rectangular contour in the complex p 4 plane having vertices in 
(— oo, oo, oo — ifi, — oo — ifi). As the poles © are defined in Minkowski space, in Euclidean space 
this means that the only poles that enter the contour have Re(£fc) < fi. Therefore 



d P4 i o ■ s dV(ip 4 



2tt v " x "*' dpi V(ip 4 ) 



_^ OO 

— = -2tt^ [Res(-i£ fc ) + Res(-i££)] 9(fi - e k ) . 



k=l 



The residue can straightforwardly be shown to be l/2iri in both cases, giving 

n(fi)=iN f N c f^J -0^8(fi~e k ). 
The theta function defines a "generalized" Fermi momentum 



(24) 



(25) 



(m«)21 



(26) 
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The quark number density can now be obtained as 

2N f N t 



3tt 2 



(27) 



fe=i 



It is actually remarkable that the density thresholds depend only on the real parts, m^. The 
imaginary parts m£ enhance the Fermi momenta (and thus the density) compared to the values 
one would get for to£ = 0. For m£ > and 

1 1 1 

< 



M 2 " (m[? : ) 2 {ml) 2 

the Fermi momentum is even larger than fi. This point will later become important. 
Integrating the expression 



fikin(0,/i) 



the thermodynamic potential can be reconstructed in a closed form 

2iY/iV c ~ 



^kin(0, fl) 



3tt 2 



fe=i 



where 



= - f^W + 5(y 2 - + 2 M 4 ] + A (a; 4 _ 6x 2 y 2 + yi) bg 



y 2 - x 2 + 2^i(p F + fi) 



x 2 + y 2 



-xy(y 2 — x 2 ) arctan 



2xyp,p F 



(x 2 — y 2 )/i 2 + 2x 2 y 2 



(28) 



It is straightforward to see that in the case of only a pair of real mass poles ±m we get the familiar 
expression for the free, massive, relativistic Fermi gas 



^kin(0,/U) 



N f N c l 



2p?pp — 5m 2 fipF + 3to 4 log 



Pf + ft 



37T 2 8 \ Dl 

We mention that at zero temperature the PL decouples, so it has no effect on the EoS. 



(29) 



III. INSTABILITIES IN A NON-LOCAL CHIRAL QUARK MODEL 

The dispersion relations Ek which enter Eq. (|14p are governed by the analytic structure of the 
quark propagator, so that further insight can be obtained only by studying the thermal behavior 
of the quark propagator, i.e., by understanding how the CCMPs respond to a change in the 
temperature or density. 

Parametrizing the analytic structure, say, from lattice studies at finite T is very demanding. 
In this case the analytic structure is also somewhat arbitrary as the quark propagator is known 
only at a finite number of points, allowing for different meromorphic forms [32j j . For the present 
purpose we will therefore study a specific model as an example case. More precisely, we consider 
a Dyson-Schwinger model with a separable gluon interaction [l2l. [l6l [l7|. In the rainbow- ladder 
approximation, these models are in fact identical to mean- field non-local NJL models, see, e.g., 
(2lL 0, Hll . They capture the important aspect of momentum dependent dressing functions in the 
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quark propagator ^ by introducing regulator functions which also ensure the convergence of loop 
integrals. Here we consider the particularly simple rank-1 case, where A = C = 1, while 

B(p 2 )=m + bf (p 2 ) , (30) 

with b being the chiral symmetry breaking parameter (mass gap) and fo{p 2 ) the regulator function. 
The latter is an input of the model. 

For this kind of separable models, it was already observed in the literature that pressure in- 
stabilities appear in certain regions of the T-/i plane. In Refs. [l3l - [l6l . l54j this was found in full 
numerical studies, and in [54| also by restricting the calculations to a finite number of CCMPs. 
The aim of the present section is to demonstrate that these instabilities are driven by the presence 
of the CCMPs, and then to study the effect of the PL. Results from a full numerical study will 
be confronted with a calculation where we restrict ourselves only to a finite number of CCMPs, 
demonstrating that the instability region is actually completely dominated by the first quartet. 



A. Analytic structure 

The analytic structure of the model was detailed in Refs. [l3l - [l5| for Gaussian and Lorentzian 
regulators. We briefly summarize their analysis for the case of the Gaussian regulator, given by 

/o(p 2 ) = e V/A » , (31) 

with a parameter Ao. The quark propagator has then an infinite number of CCMPs, as exemplified 
on the left plot in Fig. [TJ The position of the poles is controlled by the value of the gap. If the 
gap b is larger than a critical value b c given by 

b c = \ (^Jm 2 + Ik 2 - m\ e -^/™ 2 +^V^l , (32) 

all the poles are complex. For the quartet nearest to the origin, an especially interesting situation 
occurs. As the gap gets smaller, the poles travel to the real axis, where they meet in doublets at 
b = b c . If the gap is further reduced, b < b c , every doublet again splits, with one pole eventually 
going to plus (minus) infinity and the other arriving at m (—m), for 6 = 0. At that point, the real 
parts of the higher quartets go to infinity, while the imaginary parts go to zero. See Fig. Q] for the 
behavior of the first and the second quartet w.r.t. the mass gap. 

In our numerical calculations we adopt the parameters of Ref. [Hj], Ao = 0.687 GeV, m = 
0.0096A , and D = 128/Aq, successfully reproducing low-energy phenomenology. Here D is the 
strength of the non-local effective gluon interaction (i.e., four-quark non-local NJL interaction). 
For these parameters one obtains b c = 0.295 GeV, while the vacuum solution of the gap equation 
is b V ac = 0.678 GeV. Thus, the gap is overcritical in this case. 



B. Thermodynamic potential and in-medium mass gap 

The in-medium properties of the model are obtained from the thermodynamic potential ([1]), 
coupled to the PL. In mean-field approximation Q it is given by 

n(b, $, $) = n cond (b) + ft ki „(6, $, *) + §) , (33) 

see, e.g., [la ]. Here f2 C ond = Nf-gj^b 2 and U represents the mean-field PL potential, for which we 
use the familiar polynomial form found in Ref. [5l|. Other forms of this potential are in use, like 
the logarithmic for m 15211 . a strong-coupling inspired one [53j or a /i-dependent one [55j . For recent 
developments, see |56T - l58l |. 
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Re(z) b [GeV] 

FIG. 1: (Color online) Left: The curves where the real (blue, solid) and imaginary (red, dashed) parts of the 
propagator denominator T) = — z 2 + B 2 { — z 2 ) vanish for the value of the mass gap b in the vacuum. Right: 
Real parts (blue, solid) and the corresponding imaginary parts (red, dashed) of the quark propagator mass 
poles as functions of the gap parameter b. The lowest lying poles correspond to the thick lines, the next 
higher lying quartet to the thin lines. 



f2kin is provided by Eq. ([7]) augmented with the PL. This amounts to 

+°° r ^3 

fi kin (&, $, $) = -2N c N f T / 7^4tr c log [f n + B 2 (p 2 n )] , (34) 

where we now understand as diagonal matrices in color space 

Pl = P 2 +^l, w„(0 3 ,0 8 ) = oj n - + A 3 3 + \ 8 cf) 8 . (35) 

To simplify the calculations, a further restriction is imposed by setting <j)% = 0, i.e., $ = $. 
The thermal properties of the model then follow from the minimization of the thermodynamic 
potential w.r.t. b and 03 . In particular, the T and fi dependence of the dressing function B is 
solely determined by the mass gap b, as we have seen above. The explicit form of the mass gap 
equation is presented in the appendix. 

In Fig. [5] the mass gap of a system without PL is shown as a function of the temperature along 
different lines of constant n/T. The results of a full numerical solution (full lines) are compared 
with approximate ones where only a small number of poles is taken into account. The dotted 
lines indicate calculations where the system is approximated by only the lowest lying quartet, 
containing the states that become physical quark degrees of freedom when b drops below b c . In 
order to demonstrate the convergence towards the full numerical results, we also show the effect 
of additionally including the second quartet (dashed lines) . 

An observation that will be crucial later on is that, at low temperatures, a perfect agreement 
with the numerical solutions is obtained already with the first quartet. Deviations start only after 
chiral restoration, so that at higher temperatures higher quartets are needed to develop the correct 
chiral behavior. In fact, for any finite number of poles, the mass gap increases again after reaching 
a minimum, so that the correct high-temperature limit is only reached if all poles are included. 

As seen in Fig. G2 similar conclusions hold when the PL is introduced, although the deviations 
are slightly more pronounced after the chiral/dcconfincment transition, and for higher chemical 
potential also around the transition. Let us observe that at [ijT = 1 (rightmost plot in Fig. [3]), 
the system develops a first-order transition. 

In Fig. [3]we also show the results for the PL expectation value. While the overall behavior is the 
expected one, rising from $ = at low T towards $ = 1 at high T, it turns out that there is an 
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1 It .... i .... i .... i .... i .... a u.u it .... i .... i .... i .... i .... n u-u it .... i , ... i , ... i .... i ... , Tl 

0.0 0.05 0.1 0.15 0.2 0.25 0.0 0.05 0.1 0.15 0.2 0.25 0.0 0.05 0.1 0.15 0.2 0.25 

T [GeV] T [GeV] T [GeV] 



FIG. 2: (Color online) The mass gap as a function of temperature along lines of constant fi/T, for a 
system without the PL. The full line is the complete numerical calculation, while the dotted and the 
dashed lines correspond to the approximation accounting for the first and the first plus second quartet of 
poles, respectively. 




0.0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.0 0.05 0.1 0.15 0.2 0.25 0.3 0.0 0.05 0.1 0.15 0.2 0.25 0.3 

T [GeV] T [GeV] T [GeV] 



FIG. 3: (Color online) The mass gap (blue), and the PL (red) as functions of temperature. Line styles as 
in Fig. H 



intermediate regime below the cross-over temperature, where $ gets negative. Although formally 
allowed by the definition (JTTJ), which at the mean- field level boils down to $ = -^-[1 + 2 cos (/3<ps)], 
it is in sharp contrast to the standard interpretation of $ as exponential of the free energy F q 
of a static color source, $ = e~^ Fq 59]. We also note that, roughly in the same region, the gap 
parameter b rises as a function of T, a feature, which is even more pronounced in the calculation 
without PL, Fig. [2] 
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C. Pressure instabilities and instability suppression 



To calculate the EoS, the kinetic contribution to the thermodynamic potential is regularized by 
subtracting the zero-point energy of free quarks, i.e., 



"5 = - H ^/l? l0S 



—2 

p 



and the EoS is given by p(T,fi) = — f2(T, //) — tto, where is a constant chosen to achieve zero 
pressure in the vacuum. 

In Figs. |4] and [5l the pressure is displayed as a function of temperature, again along lines of con- 
stant fi/T. In the model without PL, the results are scaled by the pressure of N c Nf noninteracting 
massless quarks, 



P q sB = N c N f 



"7tt 2 


1 




180 " 


"6 





1 //^ 



12tt 2 \T 



T 4 , (37) 



whereas in the case with PL, we divide by the full Stefan-Boltzmann (SB) pressure of N c Nf 
massless quarks and iV 2 — 1 gluons, 

Psb=^b + (^c 2 -1)^ 4 - (38) 

The results for the model without PL are displayed in Fig. [4] The most striking features are 
the oscillations, which signal the thermodynamic instabilities, we have anticipated from the cosine 
terms in Eq. (|16p . They turn out to be particularly troublesome, as there are not only temperature 
regions where the pressure drops with increasing temperature, but where it gets even negative.* 
Comparing the three panels of the figure, the results seem to be rather independent of the ratio 
fi/T. We should keep in mind, however, that the pressure is scaled by the SB value, which is 
larger for larger values of fi/T '. Taking this into account, the instabilities grow with the chemical 
potential, since the Boltzmann factors are even less effective in damping the oscillating terms. This 
results in a rather large negative pressure for fi/T = 1. 

For comparison we show again the results obtained when we only take into account the lowest- 
lying poles. In agreement with our findings for the mass gap, we observe that, at low temperatures 
and more importantly, in the region of the instability, the pressure given by just the first quartet is 
an excellent approximation. The oscillations of the pressure can thus be understood quantitatively 
from the temperature dependence of b shown in Fig. [5J together with the b dependence of the 
lowest-lying poles shown in the right panel of Fig. Q] In particular, at T around 150 MeV, the 
mass gap drops below b c , so that the lowest quartet splits into two real doublets and no longer 
yields an oscillating behavior. 

At high temperatures, the full numerical result for the pressure (solid lines) approaches the 
SB limit, whereas the restriction to the first quartet (dotted) strongly overshoots this limit and 
is, thus, not a good approximation in this regime. The inclusion of the second quartet (dashed 
line) leads to some improvement but fails as well to reproduce the SB limit. This is consistent 
with Fig. [2J where the restriction to a few mass poles even qualitatively failed to reproduce the 
high-temperature behavior of the mass gap. 

Introducing the PL leads to a dramatic improvement of the EoS. As demonstrated on Fig. [5J 
the oscillations are strongly suppressed. Since the PL does not eliminate all the cosine terms 
completely (see Eq. (|21|)). residual wiggles are still present on the results for fi/T — 0,1/2, while 
at n/T = 1 also a slightly negative pressure is observed in the full numerical calculation. We also 
note that the negative values of the PL, which we have seen in Fig. [3J appear roughly in the same 



* Let us recall that we defined the vacuum pressure to be zero, so this result is clearly unphysical. 
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FIG. 4: Scaled pressure p/p q g B as a function of temperature, for a system without the PL. Line styles as 
in Fig. [3 
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FIG. 5: Scaled pressure p/psB as a function of temperature, for a system with the PL. Line styles as in 
Fig.H 



temperature region, but are more pronounced at low chemical potentials. So, to some extent, there 
seems to be a trade-off between an unphysical behavior of the pressure and an unphysical behavior 
of the Polyakov loop (when interpreted as exponential of the free energy of a static quark) . 

D. Effect of the poles at T = 

In the final part of this section, we want to discuss another consequence of the CCMPs, which 
shows up in the non-local chiral model at zero temperature. We recall that the PL decouples at 
T = and therefore has no effect in this case. 

In the left panel of Fig. [51 we show the regularized thermodynamic potential as a function 
of the gap parameter b for several values of the chemical potential ji. According to Eq. (|27p. 
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FIG. 6: (Color online) The first plot shows the behavior of the regularized thermodynamic potential at 
T — as a function of the mass gap. The curves correspond to different chemical potentials: /i = (solid), 
H = 0.45 GeV (dashed), /i = 0.52 GeV (dash-dotted), and /i = 0.58 GeV (dotted), where chiral symmetry 
is broken again. This is made transparent on the second plot, where the mass gap is displayed as a function 
of the chemical potential. On the third plot the resulting EoS is shown (full line), in comparison with the 
EoS for a gas of free gas of quarks (dashed line). 



the threshold chemical potential for non-zero quark number density is controlled by the lowest 
value of mj? . Hence, at finite chemical potential, the thermodynamic potential has to stay at 
its vacuum value in those regions, where the lowest threshold, m^, is bigger than ii. This is 
similar to the local NJL model, where at T = the thermodynamic potential as a function of 
the constituent quark mass M stays at its vacuum value for M > fi.> In the non-local model the 
essential difference is that the threshold is a non-monotonic function of the gap parameter 
b. As we have seen in Fig. [TJ it rises from = m at b = to a maximum at b — b c , where 

m i (b c ) = | (rn + \J m 2 + 2Aq^ = 0.489 GeV for our parameters. Above this point, mf slowly 

decreases, reaching arbitrarily small values at large b. As a consequence, for < fi < m*(6 c ), 
the thermodynamic potential coincides with the vacuum curve only in a finite interval around 6 C , 
bounded by the condition mf (b) — fx. For the mass gap outside this range there is a finite density 
of quarks and the thermodynamic potential is below the vacuum one. An example for this case is 
given by the dashed line in the left panel of Fig. [SJ 

For ix > mf(6 c ) the potential is reduced everywhere (dash-dotted line). Eventually this leads 
to a pathological result: Whereas at intermediate chemical potentials a global minimum emerges 
near 6 = 0, leading to (approximate) chiral-symmetry restoration, at sufficiently high chemical 
potential, the nontrivial minimum at large b becomes the global one, meaning that chiral symmetry 
is broken again (dotted line) . The corresponding behavior of the gap parameter which minimizes 
the thermodynamic potential is shown in the middle of Fig. [5] This obviously unphysical result is a 
consequence of two facts: First, for fi > mf (b c ) there are two CCMPs, i.e., twice as many "degrees 
of freedom" which contribute to the pressure in the non-trivial minimum. Second, according to 
Eq. ([26]) . the Fermi momentum is lowered by the real part but enhanced by the imaginary part 
of rofc, and eventually becomes even larger than [i. These effects lead to a further enhancement 
of the density and thus the thermodynamic potential decreases faster in the non-trivial minimum 



t More general, this is a consequence of the so-called "Silver Blaze problem" |60H . meaning that at T = the grand 
partition function must not change if fi is below the smallest excitation threshold. 
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than in the trivial one. This is also underlined by the right plot, where the pressure is displayed 
as a function of the chemical potential: In the second chirally broken phase, the immense rise in 
the density causes the EoS even to overshoot the pressure of a free quark gas. 

The second breaking of chiral symmetry at high chemical potential is reminiscent to what we 
have found at high temperature, when we restricted ourselves to a finite number of poles, cf. Figs. [5] 
and [3] In fact, at T = 0, the chemical potential imposes a sharp cutoff on the thermal part of the 
thermodynamic potential, and in this sense, there is always only a finite number of poles which 
contribute. For instance, the highest value of \i considered considered in the left plot of Fig. [5] is 
still smaller than for the shown range of b, and therefore the thermodynamics is determined 
completely by the first quartet. 

IV. SUMMARY AND CONCLUSIONS 

A realistic continuum effective theory of strongly interacting quarks and gluons should incorpo- 
rate the two most striking features of low-energy QCD: chiral symmetry breaking and confinement. 
The first one represents a physically familiar concept, and it is rather easily modeled. The latter 
is still lacking a proper explanation, but can nevertheless be realized by various confining criteria 
like, e.g., positivity violation. 

The simplest possible model of chiral quark dynamics is the Nambu-Jona-Lasinio model, exhibit- 
ing real mass poles, a feature shared with its nonlocal, but instantaneous generalizations. Although 
very successful, these models do not incorporate confinement. Models in which the quark propa- 
gator is an entire function represent the opposite situation that quasiparticle poles for quarks (and 
gluons) are absent, but therefore do not allow for deconfinement. In between, there is a large class 
of models that have poles in the complex plane, which, due to the their manifest covariance, come 
in complex conjugate pairs. 

In this paper we have reported the observation that the thermodynamical state in such a class 
of models is unstable, due to the possibility of oscillating, or even negative, pressure. Further con- 
sequences can also be entropy decrease with the temperature or negative heat capacity, violating 
the standard stability criteria for thermodynamical equilibrium. This underlines that quark con- 
finement is a multifaceted phenomenon that cannot exclusively be modeled as strong momentum 
dependence of dynamical mass and wave function renormalization. In this context we have inves- 
tigated to what extent a suppression of these unphysical instabilities can be achieved by coupling 
the system to the Polyakov loop. 

A simple quark model which reproduces the CCMP form of the quark propagator dynamically is 
a DSE model with a covariant separable interaction. As an example, we solved such a model for a 
Gaussian formfactor ansatz at finite temperature and chemical potential in mean-field approxima- 
tion in order to demonstrate the possible effects of CCMPs on the thermodynamics. The results 
show that CCMPs are indeed the mechanism responsible for the instability. When coupling the 
quark propagator to the PL, we find that the pressure instabilities are strongly reduced. Unfortu- 
nately, more or less in the same region of the phase diagram, the PL itself becomes negative, which 
is in sharp contrast to its standard interpretation as exponential of the free energy of a static color 
source. Thus, although there is no exact one-to-one correspondence, it seems that one unphysical 
effect could only be suppressed to the expense of a new one. 

At zero temperature and finite chemical potential the PL is irrelevant. For very high chemical po- 
tentials (~ 500 MeV in our case) CCMPs produce yet another unexpected and probably unphysical 
behavior: the chiral symmetry gets broken again. * We expect, that the same behavior is persistent 
in a more realistic setup, i.e., when one includes the effects of wave function renormalization as, 
e.g., in 



f At very high chemical potential chiral symmetry is expected to be broken again in a color-flavor locked phase .66] , 
but this is a completely different mechanism. 
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From a wider perspective, vacuum instabilities in a mean- field (or "classical") description of 
QCD are not uncommon; a prominent example being the tachyonic mode observed in the Savvidy 
vacuum [62|,[63|]. It is interesting to note that also in this case the PL acts as a "stabilizer", i.e., by 
suppressing the original tachyonic modes |64j. 65] . Unfortunately, this program is also not entirely 
successful as new unstable modes arise [64 165|. 

Apparently, the lesson to be learnt is the following: in order to cure the thermodynamic insta- 
bility problems in the low-temperature and low-density domain of effective "confining" models of 
quark (and /or gluon) matter properly one has to go beyond the mean-field level of description, 
however clever it may be designed. Only the explicit inclusion of the physical degrees of freedom 
in that domain, the hadrons as color neutral bound states of quarks and gluons, will provide the 
non- negative and monotonously rising pressure in the confinement domain of low-energy QCD. As 
first promising works in this direction we refer to [ljl [U [23|, [26| and suggest to develop the CCMP 
propagator class models beyond mean field. 
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Appendix A: Gap equation in medium 



This appendix presents the quark gap equation in the mean-field approximation, with the Mat- 
subara summation being analytically performed. 

Minimizing the mean-field thermodynamic potential (|33[) . the quark gap at T, fi > is obtained 
to be 



16-Dp 
9 



T 



+ 00 

E 



d 3 p 



tr c 



B(P 2 n )fo(pl) 



(Al) 



where p n is to be understood as a diagonal color matrix, see (1551 . The sum over Matsubara 
frequencies is evaluated using the standard technique of rewriting it as the sum over residues of 
a contour integral in the complex energy plane over the analytically continued integrand function 
folded with the function (|18p having simple poles at the PL shifted Matsubara frequencies. The 
Matsubara summation is thus converted into three closed contour integrals, similar to what was 
performed in Sec. [TXJ These are calculated by the residue theorem, giving the result 



b = 



16 Ai 



Nr. 



d A P B(p 2 )f a (p 2 ) 
(2tt) 4 p 2 + B 2 (p 2 ) 



k=l 



d 3 p 



Re{Res(£ fc )tr c [n + (£ k ) + n-(8 k )]} 



(A2) 



The first term is recognized as the vacuum gap equation. Residues can easily be deduced to be 



Res(f fc ) = 



j?(p 2 ,-g 2 )/ (p 2 ,-£ 2 ) 
V'{£ k ) 



(A3) 
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with T> given by (fTU|) when A = C = 1. An expression for a color trace of the occupation numbers 

3$ e -/3(£fcTM) _|_ g$g-2^(£fcTM) _|_ 3 e ~3/3(£ k Tv) 



tr c [n±(£ fe ) 



1 + 3$e^^( £fc= FP) 4. 3$ e -2/9(£fcTM) 4. e -3/3(£fcTM) 



(A4) 



completes the calculation. 

At zero temperature the Matsubara sum is converted to an integral which is performed in a 
similar fashion as the integral for the quark number density, see Eq. (f24f. The result reads 



d 4 P 5(p 2 )/o(p 2 ) 
(27r) 4 p 2 +£ 2 (p 2 ) 



2E 

k=l 



d 3 p 
(2^)3 



Re{Res(£ k )}0(ii-e k ) 



(A5) 



with Res(£/c) given by Eq. (|A3 
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